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Crosslinked semi-flexible and flexible filaments that are actively deformed by molecular motors 
occur in various natural settings, such as the ordered eukaryotic flagellum, and the disordered 
cytoskeleton. The deformation of these composite systems is driven by active motor forces and 
resisted by passive filament elasticity, and structural constraints due to permanent cross-links. Using 
a mean field theory for a one-dimensional ordered system, we show that the combination of motor 
activity and finite filament extensibility yields a characteristic persistence length scale over which 
active strain decays. This decay length is set by the ability of motors to respond to combination of 
the weak extensional elasticity, passive shear resistance and the viscoelastic properties of the motor 
assembly, and generalizes the notion of persistence in purely thermal filaments to active systems. 


Semi-flexible and flexible filaments interacting with ac¬ 
tive molecular motors arise naturally in the cell as well as 
in in-vitro motor assays [THIS]. The overall coordinated 
macroscopic mechanical response of such a composite as¬ 
sembly depends on the coupling between the active mo¬ 
tors and the passive properties of the aggregate compos¬ 
ite. A structurally ordered example of such a system is 
the eukaryotic flagellum comprised of relatively stiff fila¬ 
ments (microtubules), motors (dynein) and passive elas¬ 
tic elements (nexins) which together oscillate with well 
defined wavelengths and frequencies [HI [T4HT6] . To un¬ 
derstand how coordinated motion can arise in such a sys¬ 
tem, a variety of models of varying degree of complexity 
[9HI3] have been proposed; those show that the wave¬ 
length of active beating is correlated with increasing flag¬ 
ellum length. However, experimental evidence strongly 
indicates that this relationship is altered for very long 
lengths nzmi] and that the wavelengths attained in na¬ 
ture are self-limiting m even as the flagella themselves 
range from tens of microns to nearly a centimeter. This 
strongly suggests that over large lengths, mechanical in¬ 
formation transmission degrades substantially. 

Indeed, in a passive context HU, it has been demon¬ 
strated that when bundles of filaments are forced to bend, 
the shearing forces between them are mediated by ex¬ 
tensibility, leading to a characteristic scale over which 
mechanical information is transmitted. In an active sys¬ 
tem such as an array of soft filaments interacting with 
each other via motors, degradation of mechanical signal 
transmission occurs because the weak passive extensibil¬ 
ity of the filaments limits the range over which motor- 
co-ordination can occur, while restrictions on the kinet¬ 
ics of motor activity set by ATP hydrolysis rates couple 
back to weak extensibility and alter the interaction be¬ 
tween motor groups significantly. We illustrate this in 
Fig. 1(a) showing two distinct motor patches labeled I 
and II that cross-link a pair of thin filaments that may 
bend and shear. For inextensible filaments, mechanical 
activity by group (I) leads to shear (sliding) in this patch 
which is transmitted by attached motors in (II) over ar¬ 
bitrarily large inter-aggregate distances. For extensible 
filaments, there is an interplay between shear/slide and 



FIG. 1. (a) An active, motor-filament composite illustrating 
two distinct motor aggregates, groups of motors (I) and (II). 
Of interest is the distance between two motor groups beyond 
which they act as if isolated, (b) A sketch of the coarse¬ 
grained active gel (top) corresponding to the mean field ap¬ 
proximation of the microscopic filament / motor system (bot¬ 
tom) illustrating the motors (blue) and cross-linking passive 
linkers (red). 


filament elongation along its contour and this leads to a 
decrease in the ability of group (I) to mechanically link 
with group (II). In this letter, we propose a minimal, 
mean-field theory to study the decay of deformations gen¬ 
erated by localized steady or oscillatory strains in active, 
weakly noisy motor-filament assemblies. For simplicity 
we focus on purely extensional systems ignoring bend¬ 
ing modes of deformation and restrict ourselves to a one 
dimensional setting. 

The model system we analyze is illustrated in Fig. 
1(b). A weakly extensible filament of length £, thick¬ 
ness w Young’s modulus E and lateral extent b is 
animated by active motor proteins that are spaced uni¬ 
formly with fixed areal density These motors are at¬ 
tached to an underlying rigid, straight substrate. Passive 
permanent cross-linkers connect the filament to the sub¬ 
strate maintaining a gap of size D {w < D <^ £) between 
the filament and the substrate. These cross-links, spaced 
uniformly with areal density also act as linearly elas¬ 
tic springs of stiffness Active motors are modeled as 
linear elastic springs of stiffness km^ that attach and de¬ 
tach from the elastic filament. The inter-link spacing and 
the inter-motor spacing is assumed to be much smaller 
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than the filament length so that a continuum model is 
appropriate. Thus when b ^ max(re,-^), the assembly 
acts as an effective one dimensional extensible composite 
filament, with the passive permanent cross-linkers that 
connect the filament to the substrate providing shear 
resistance. The effective shear modulus resulting from 
these passive springs is obtained by considering the rel¬ 
ative displacement between the top sheet and bottom 
sheets due to a shearing stress acting on an area b£; this 
yields G* ~ and CwjD ^ wkf^pf^/2 - consis¬ 

tent with analysis of a passive cross-linked railway track 
model [To]. 

Due to the motors, the filament is displaced horizon¬ 
tally as it is stretches. Let the displacement of a mate¬ 
rial point on the filament relative to its initial position 
in the undeformed state be [/(s,t), s being the mate¬ 
rial/reference coordinate, and t time. For an inertialess 
filament, the sum of forces on the filament arising from 
filament extensibility, springs and motors must vanish. If 
Fm is the effective active force density per attached mo¬ 
tor and pm the (total) motor density, local stress balance 
then yields 

ds{Kpa.sdsU) — GpasG + pmFm = 0- (1) 

Here LCpas ^ Fw and Gpas ^ Pn^n are the passive 
stretching and passive shear moduli (per unit width) 
for the composite filament. In the absence of activity, 
Fm = 0 SO that the solution to (1) yields an exponen¬ 
tial profile for the displacement with a relaxation length 
scale Ie = (^pas/Gpas )2 = (F^ic/pn^n) 2 , so that as the 
passive shear resistance Pn^n increases, the persistence 
length decreases. 

When motors are present, Fm 7 ^ 0. At the level of a 
single motor, we assume that motors attach to the elastic 
filament with an initial extension dm (that may vanish) 
and are then stretched due to two effects - a sliding ve¬ 
locity relative to the filament as well as the displacement 
of the filament due to the cumulative effect of other mo¬ 
tors. Both of these lead to an increasing motor extension 
y shown in Fig. 1 (b). Assuming a simple linear relation 
for the force-velocity relation for an individual motor, we 
may then write dty = dtU Vo{'\- — kmy/Fg) where the 
first term is the velocity of the filament and the second 
term is the velocity of the head relative to the track. Here 
the stall force Fg and the load-free velocity Vq are motor 
dependent parameters. 

With this given motor-filament interaction, the aver¬ 
age number of attached motors relates attached and de¬ 
tached probability densities to their corresponding fiuxes 
via microscopic transition rates. In the mean-field limit, 
the density of motors pm is large enough that fluctua¬ 
tions in the time average of the density of attached mo¬ 
tors pa = pmN{F)^ where N{t) is the fraction of attached 
motors are small compared to the mean value. Ensemble 
averaging the terms in the population balance equation 


over y (see SI-§I) yields 

dtN = (2) 

where the mean attachment rate is assumed to 
be constant, and the mean detachment rate cjoff = 
is modified by motor extension, since the 
characteristic energy scale for detachment is £ = 

km^m/^BT and depends on a typical length scale Sm 
charactering the extension at which motors detach. Con¬ 
sistent with previous studies [lain], we use two forms: 
F{£,Y) = exp(f|F|) and F{£,Y) = cosh(fT). The 
evolution of the mean motor extension Y = {y)/Sm is 
obtained by write down the microscopic description for 
dty for a single motor and then average over y. Relating 
the average stretching rate {dty) to dt{y) yields 

dtY = ^+(^2 - i")+<„(^3 - i") (^r) 

(3) 

The first term on the right is the passive convection of 
the motor head as the filament stretches m and sec¬ 
ond term arises as the motor head walks relative to the 
moving filament. The last term on the right corresponds 
to the rate at which the mean strain changes due to the 
kinetics of motor attachment and arises from the differ¬ 
ence in extension of attaching and detaching motors. The 
parameters Ai = Vokm/(^off^s), ^2 = Fs/{kmSm) and 
As = dml^m relate the microscopic motor parameters to 
the population level dynamics. Finally, given the frac¬ 
tion of attached motors N and their mean extension F, 
the effective active force per motor in 0 is given by 

Fm = kmSmNY. (4) 

Equations (l)-(4) specify the evolution of the filament 
displacement U{sA)^ the motor density N{sA) and the 
motor stretch Y{sA) in terms of the dimensionless pa¬ 
rameters Ai^A 2 ,As and T = While T is re¬ 

lated to the duty ratio of the motors, Ai and As deter¬ 
mine the mechanism that yields the average extension of 
attached motors. Specifically, when ^ 3 ^ 1 , the role of 
the pre-strain dm is negligible m while when Ai 1 , 
pre-strain dominates [ 12 ]. Our model complements and 
differs from previous attempts in a few important ways. 
Eirstly, the detachment rate depends on the mo¬ 
tor extension [laiiH] and not the rate of extension m, 
consistent with experiments that show that bond failure 
is more naturally dependent on the extension and only 
weakly on the rate of extension [19]. Secondly, the non¬ 
linear coupling between passive and active deformations 
in 0 distinguishes our model from previous studies of 
motor mediated bending of filaments [18]. Einally, we as¬ 
sume that the filament is (weakly) extensible, bolstered 
by experimental evidence in passive systems lauB, and 
unlike all previous studies of active ordered filaments. 

To understand the solutions of (l)-(4) and their sta¬ 
bility, we expand the motor displacement F, the fila¬ 
ment extension I/, and motor number density N in pow- 
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ers of e ^ 1, the deviation from steady state, using the 
ansatz {U,Y,N) = ([/q, Iq, A^o) + e Real[e^^(^, F, TV)], 
with the first term corresponding to the steady, static 
state (that arises due to motor activity) while the second 
corresponds to the time dependent perturbations from 
the steady state - these could be due to intrinsic or im¬ 
posed perturbations. Here €Uimp/Uo{i) 1 and a is a 
complex growth rate. In the absence of external forcing, 
this growth rate determines the stability of the base state 
(I/q, Yq^Nq). Substituting in 0-(|^ we find at O(e^) that 
Uq.Nq.Yo satisfy 

~ Gpasb^O + Gg,ctC = 0, (5) 

TVo = (l + ^J’o)-', ( 6 ) 

1^0 = (.41^2 + A^Ta){Ai + J-o)-' (7) 

where Gact = pm^mi the active analogue of passive shear 
modulus and C = (^^AToYo = • Here, 

No and Fq are independent of position, s. The steady 
extension at the free end Uo{t) results purely from steady 
motor activity when e = 0 and is therefore obtained by 
imposing the constraint dsUo{i) = 0. Equations (5)-(7) 
determine the extensional (strain) field resulting due to 
activity. At 0(e), we find 


F 


T^pas OssU GpasH -|- Gact 5m{NoY + YoN) = 0, (8) 

N{a + uj°A^ + '^Ao))+co°r,^NoKY = 0, (9) 

^ + (^3 - Yo)^ - = 0 - ( 10 ) 

^on / -‘■^0 (Jm^on 


Equations (8)-(10) may be interpreted in two ways. In 
the absence of imposed oscillations, they determine the 
linear stability of the equations (5)-(7) to small pertur¬ 
bations. Alternately (8)-(10) can be interpreted as the 
linearised response to a small amplitude, externally im¬ 
posed oscillation in which case in the limit t ^ oo limit 
(provided Real [a] < 0), (G, TV, F) is slaved to the forcing 
frequency. 

Before we approach the question of persistence of ac¬ 
tivity in an active filament, we ask if spontaneous strain 
fields can be generated by a small, localized patch of mo¬ 
tors of length is ^ (A"pas/Gpas) 2. The dynamics of this 
patch may be mapped to that of a rigid segment working 
against an effective spring with passive as well as active 
components - here, we lump these contributions together 
using an effective spring constant Kg ^ Gpas (SI-§Ha). 
The displacement of this fragment follows 


-iV,G + Gact^^TVF = 0. (11) 


When Ks > 0 (for the degenerate case = 0, see SI- 
§Hb) and global filament translation is prevented, equa¬ 
tions (2), (3) and (11) admit two classes of solutions - a 
constant extension of the rigid filament with dtU = 0, or 
a time-periodic, oscillatory extension. In the first case, 
the steady density and motor strain are given by equa¬ 
tions (6) and (7) with the steady extension GactG/iVg 
where G = (^^TVoFq. This static state is however lin¬ 
early unstable to oscillatory states and the onset of these 



FIG. 2 . Oscillatory solutions for A 3 = 0, As > 0 and 
Apas = 0 obtained by analysis of 0 - (a) Neutral stabil¬ 
ity curves T = vs. number of motors (proportional 

to motor density) are shown for S = 4.944, Yo = 5 x 
and ccon = 0.25 s.“^ Curves are - (1) = cosh(^F) and 

(2) T — exp(^|F|). For fixed T, there is a critical num¬ 
ber of motors required for oscillations to manifest (as illus¬ 
trated for T = 0.6). There is a critical value of Kg I km for 
fixed motor density below which oscillations will not be seen, 
(b) Scaled frequency of oscillations, ujcI^on cts a function of 
T. The solid line is the exact result while the dashed line is 
the adiabatic approximation Ft TFJ^. (Inset) The model 
spring-filament-motor homogeneous system. 


oscillatory states can be readily obtained from (9), (10) 
and the linearised version of (11). Analyzing the condi¬ 
tions under which Real [a] = 0 yields the neutral sta¬ 
bility curve that demarcates the static state from sta¬ 
ble, oscillatory states, emanating via supercritical Hopf- 
Poincare bifurcations; the necessary criterion for this be¬ 
ing ^(Fq) = > 0. The mechanism behind the 

oscillations is the coupling in (11) between the fila¬ 
ment motion and attached motors that yields an effec¬ 
tive compliance comprised of both in-phase and an out- 
of-phase term; the oscillatory solutions obtained are a 
result of positive feedback, when either or both of these 
terms turn negative, similar to that seen in flagellar m 
and spindle oscillations [12]. In FigJ^(a), we show neu¬ 
tral stability curves (SI-§HC) separating oscillatory from 
non-oscillatory steady regimes. For fixed motors den¬ 
sity (fixed motor number) there is a range of T over 
which oscillations can be sustained. For fixed T there ex¬ 
ists a critical motor density, related to Kg/{pm^m) (units 
m“^), above which the oscillations appear. When stably 
oscillating, the power due to activity balances dissipation 
due to motor friction - this yields an effective viscosity 
^act ^ pmNokm{d{y)/dUt)ut=o where pmNo is the to¬ 
tal number of attached motors. As T increases, both 
{d{y)/dUt)ut=o and Nq typically decrease leading to an 
increase in the emergent frequency uuc - c.f. Fig. ib)- 
here, Ut = dtU. 

Having identified the two stable localized dynamical 
states that can exist for a short fragment, we next com¬ 
bine this local analysis with a global description of the fil¬ 
ament extension to obtain decay lengths for active assem¬ 
blies. We first address decay lengths of steady extensions. 
For boundary conditions, we assume that A(0,t) = 0 
and U{i,t) = Uo{i) + eAiReal[e*^^], i.e. one end is fixed 
while the other is forced harmonically with amplitude 
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eUi. We focus on homeostatic responses and analyze the 
frequency-locked response of the active asssembly to the 
imposed frequency at the free end. 

At O(e^), we recover equation (5) which when solved 
with boundary conditions t/o( 0 ) = 0 and dsUo{i) = 0 , as 
is appropriate for a free end, yields 

Uo = (l-a - (1 - a) (12) 

^pas 

where a = ^ decay length 

for steady extensions depends only on the passive elastic 
compliance because active motors sense only strain rates 
and not the strain. An elastically compliant inter-motor 
spacing or inter-motor interactions will alter the decay 
length from that given by (12). In the degenerate limit 
Pn = O 5 Uo{s) = (GactC'/iCpas) s{2£ — s) - there is thus no 
decay length scale for the steady extension in this case. 
Attached motors in rigor do sense absolute displacements 
and strain rates. The decay length is then modified to 

(SI-§IIIA) (Arpas/[GactA'o + Gpas])^ where A^o is the pre¬ 
existing attached fraction. 

Next, we examine equations (8)-(10) to discern the de¬ 
cay of oscillatory strains. The end at s = is forced har¬ 
monically at frequency uj and so as t ^ oc, the locked-in 
response is obtained by setting a = iuj in equations (9) 
and (10). Combining ( 8 )-(10) then yields 

^pas dssU — GpasG — GactXG = 0 (13) 

where the complex function y = —SmiNoY - 1 - YqN)/U = 
1ZY\I characterizes the linear viscoelastic response of the 
motors (SI-§IIIB,C). Since the stationary solutions, Tq 
and A^o do not involve the position dependent Go, y does 
not depend on s. We seek solutions to (13) of the form 
U - with A being a complex wavenumber. Non¬ 
trivial solutions of (13) exist only when (Gpas + Gact^ ~ 
A^pas/^^)^ + — d which on simplifying yields the 

two roots = [{Gp!,s/Kp^s) + (Gact/-f^pas) ± 
Defining £e = \/Kpiis/G^ and an activity parameter 
P = Gact/Gpas, we find an effective persistence length 
scale, Ae = Real [A] to be given by 

\ — r 

XE = V2iE (v^rTFW^TW + i + /37e) ^ (14) 

Here /3 controls the importance of active to passive shear 
stiffening. When /3 = 0, we recover the result for a pas¬ 
sive assembly, Ae = In the absence of passive linkers, 
the persistence length is controlled solely by Gact/Apas- 
The frequency dependence of Ae, arises via the real and 
imaginary parts of y, these being associated with motor 
properties such as stall force free velocity vq and the 
duty ratio via T. When activity is weak, /3 ^ 1, an 
expansion of equation (14) in powers of p yields Ae ~ 
£e (1 - Ipn - - sn^)) -f For F > 0, at 

0{/3) the motor activity decreases Ae with shear stiff¬ 
ening driving the stretching of the filament. The 0(/3^) 


term can either enhance this effect or lessen it depending 
on the sign of the term (X^ — 3F^). For strong activity, 

when ^ > 1, we get Ae « ^j2{Kp^/G^ct){T^ + 1x1)“^ 
with actively driven shear stiffening dominating. 


The role of motor dynamics can be probed by examin¬ 
ing the low and high frequency limits of (14). In the low 
frequency limit, motors undergo multiple mechanochem- 
ical cycles before they feel the imposed oscillation. In the 
high frequency limit, attached motors encounter several 
cycles of imposed oscillations before detaching. Activity 
will then increase the persistence length from the lower 
bound of a tightly coupled system. These trends are 
transparent when we examine the limit = 0. In 
the low frequency limit, X ^ — 7 + yr) 

and IZ ^ (co’/co’off)^(l — 7 + yr^) with 7 = 
and r = <^off/(<^off + ^on) both dependent on Tq 
but not on s. Substituting these in (14) yields Ae ~ 

(1 - |(wVw 7 ){/ 32 (l -7 + 7r)2 + 4^ -7 + 7 t2)}), 

which clearly indicates that activity, through 
decreases the persistence length. In the limit of 
high cj, we find IZ 1 - (woff/w)^(l - 7 + 
and X (co’off/co’)(l — 7 + l/^) providing 

Ae ~ {£e /Vl + ^)(I + |(<^off/^)^{4/3(l — 7 + 

y/r^) — 1^(1 — 7 + 7 /^)^})- In fhe limit cj ^ oc, 
Ae ^ £eI a/I + P = \/Apas/(Gpas + Gact)- 


For fixed motor properties, the decay length sets a fi¬ 
nite range of correlated activity that might be relevant 
for naturally ordered active matter such as eukaryotic 
flagella [20]. Modeling the flagellum as a composite fila¬ 
ment driven by uniformly spaced dynein motors, we take 
re, D and E as the radius, spacing and Young’s modulus 
of microtubules and /cn and pN as the nexin stiffness (per 
length) and nexin density. From the existing literature, 
we find w ^ 20 nm, F ^ 40 nm [ 21 ], E ^ 1.2 GPa [9l|2T], 
/cn ^ 16 —100 pN pm”^, km = 10“^ N/m [22], pm ^ 10^ 
m“^ and wp^ ^ 10 ^ — 10 ^ m“^ [T71[22]. We find that the 
frequency dependent decay length £e ^ [10 — 500] pm, 
where the lower limit corresponds to all dyneins and nex- 
ins attached, and the upper bound corresponds to only 
nexins being attached - in the range of physically rele¬ 
vant length scales. For segments of the flagella separated 
by distances less than the decay length, one could have 
frequency entrainment with multiple oscillatory regions 
eventually synchronising. 


Our theory accounts for extensional deformations in a 
model ordered active composite consisting of elastic ele¬ 
ments that can be stretched and sheared by dynamically 
attaching motors. Extensibility, even if negligible locally, 
affects the scale over which deformations persist owing 
to the competition between extensibility and shear. For 
strain dependent motor kinetics, the decay of steady ex¬ 
tensions is controlled by shear stiffening due to passive 
linkers. Localized oscillatory strains lead to global de¬ 
formations that decay over longer length scales that in¬ 
volve passive and frequency dependent active properties 
of the assembly. Given that it is now possible to cre¬ 
ate partially ordered mixtures of motors and cytoskeletal 
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filaments [23l[24] that show spontaneous activity, our cal- culations can be tested experimentally. 
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SI: Persistence of strain in mot or-filament assemblies 


I. COARSE-GRAINING YIELDS MEAN FIELD MODEL 


Of the total motors per area at position s along the filament let a fraction N{t) be attached so that the density 

of detached motors is pd{sfl) = Pm(l ~ Each motor experiences an internal spring extension y which may 

be considered the distance between the head and the base or as an intrinsic co-ordinate of the underlying mechano- 
chemical cycle. We choose ^ = 0 as our base state corresponding to an unstrained motor. 

In the mean-field limit, the number of motors (and the density) is large (high) enough that fluctuations in the time 
average of the density of attached motors pa{sfl) = p^N^sfl) are small compared to the mean value. The motor 
kinetics may be described by a set of population balances relating the attached and detached probability densities 
Va{s,yfl) and Vd{s,yfl) to attachment and detachment fluxes Ja and Jd via the transition rates z^on(^) and z^off (^)5 

dtVa + dyJa = yonPd “ (15) 

dtVd + dyJd = l^oS'Pa - T^onVd (16) 

consistent with previous studies (see [H [12] and references therein). For simplicity, the attachment frequency z/on 
is chosen to be constant while the detachment rate is allowed to vary with the motor extension. Thus we write 
z^off = z^off(^ 5 ^m): being a measure of the motor extension at which the detachment rate becomes significant. 

Focusing on the noise-less limit, we make the assumption that thermal noise as described by Brownian effects of the 
diffusion of both the motor base (this is assumed to be attached permanently) and the motor head, is negligible, 
where, in the absence of motor diffusion, Ja = {dty)Va- 

Consider now the flux of attached and detached motors. One scenario is that the detached motors could rapidly 
relax to their equilibrium internal state with a viscous relaxation time. Alternately the detached motors could relax 
to a delta function Vd = pd{t)S{y) with the change occurring instantaneously. When the distribution in extension of 
attached motors is sharply peaked about the typical (average) length, transients to this distribution occur over times 
very small compared to the averaged macroscopic time scale, t we can write Va = pa{t)S{y — (y)). Averaged equations 
(moments) may be obtained from (SI-1) and (SI-2) by using the definition J XVady = {X) Npm- 

Integrating either (SI-1) or (SI-2) over all configurations y yields 

dtPa = - /5a«„ + (17) 

In deriving these equations, we have used the fact that averaging the transition rates yield the mean-field (effective) 
attachment rate (chosen to be a constant), and detachment rate, (z/off) = cjoff that is modified by the mean motor 
extension. 
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The extension of attached motors y for a location s is peaked about the mean value {y){s,t) with small deviations 
from this mean. Ignoring these deviations as done in previous studies ISIig, we choose to depend only on the 

scaled mean motor extension Y = {y)/Sm so that cjoff is just dependent on (T, ^m)- Thus we may write, 

(^off y) ~ ~ (t8) 

and thus c^off = (^ 5 ^) where S = km^‘^/kBT is a measure of the energy needed to stretch the motors to their 

maximal extent. For subsequent analysis, we define 


^o = ^(^o), 


- §(y = ^ 0 ). 


Note that since pa{s,t) = pmN{s^t)^ equation (SI-3) can be recast as 


ai7V = <„-iVK„ + u;off). 


(19) 


( 20 ) 


The final equation is for the evolution for the mean strain, {y) =YSm- This ensemble averaged quantity changes 
due to three effects - the pre-strain (pre-extension) value of attaching motors, the fact that motors detach at finite 
strain and the evolution of the strain whilst the motor is attached. The evolution of extension in attached motors is 
related to the correspondence between the change of the mean stretch, dt{y), to the mean rate of stretch, (dty) and 
necessitates consideration of the local interaction between the motor heads and the filament tracks on which they 
move exemplified by an active velocity-force relationship and a passive force that is a measure of the elasticity of a 
fixed cross-bridge under conditions of rigor. 

We consider a simple physical picture for how motors interact with the filament. We assume motors attach in a 
pre-strained state with initial extension at attachment dm and then once attached walk relative to the filament with 
velocity Vm as they interact with the tilted potential that characterises the filament-head interaction. At the same 
time since the base is fixed, the increasing spring extension y causes a force that ultimately at some point detaches 
the motor. 

As the motor head moves relative to the potential field generated by the filament, it obeys a velocity-force rela¬ 
tionship prescribed by microscopic details of the motor head-track interaction. One may then think of the velocity 
relative to the track as being caused by the hopping of the head from one potential local energy well (binding site) 
to another without detaching from the track. The simplest specification for the velocity of the motor head relative to 
the track after it is attached is given by a linear force-velocity relationship 

Vm = dty = dtU + (^1 - ( 21 ) 

where the first term dtU is the velocity of the filament (that convects the attached motor). The second term is the 
velocity of the head relative to the track and in this term Fg is the stall force at which motion stops, Vq is the velocity 
at zero force and km^mX is the spring force. Taking the ensemble average of (SI-7) yields 

/ h y(d \ 

{dty)=dtUFVo[l-^^y (22) 


Finally, we put together these aspects of the motor extension. Multiplying (SI-2) by y and then integrating over 
we can write 


dt{y) = -{vofi y) + dtU + Vo[l- 


kmY^d 


{iyof!){y) + 0 J°n{dm - (y)) 


l-N 

N 


Since (z/off y) ~ we get the dimensional version 

'l-N 


dt{y) = dtU -Uj^^{y) 


N 




^Off 


1 - 


km {y) 


+ drn 


l-N 

N 


(23) 


(24) 


Note that when dm = 0, we recover the model proposed by Grill et. al [12] wherein motors attach at zero extension 
and then walk along the filament with a strain-dependent speed. When Vq = 0, we recover a different model [6] wherein 
motors follow a two-state cross-bridge scheme and attach in a forward leaning position with an motor extension dm- 
Equation (SI-10) is an approximation which allows one to interpolate between these two limiting cases. 
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FIG. 3. Stationary branches and stability of steadily translating solutions of a fragment of the hlament when Kg = 0 and 
Ai = 0. Here, (utcOo^/Sm) = Z is the continuation parameter. Parameters values are, scaled motor extension dml^m — 1, ratio 
of kinetic, c^on/^off = 1- The function T — exp(^|y|) - solutions are shown for (a) ^ = 1 and (b) ^ = 2. Circles are limit 
points where stability changes from stable (solid line) to unstable (dashed line) regimes. Increasing the value of E opens up the 
unstable part of the branch. 


Now using Y = {y) and defining 

s) , ^2 = Fsikrn^m) : ^3 = 

we can write (SI-10) as 


dfY = 


dtU 


dm. \ JY 


l-N 




Ai^{A2-Y)yA3 


l-N 

N 


(25) 


(26) 


To find steady homogeneous solutions, we set = dfY = dfN = 0 in (SI-6) and (SI-12). This gives us 

1 


No = 


and the implicit, non-linear equation 


^0 = 


l + TJ’o 

^ 1 - 4-2 + A 3 J-Q 


Ai Y J~o 

Note that when ^ = 0, then Tq = A 3 . For later purposes we define a constant C = (5^A^oTo, so that 

1 ^ 1^2 + A3J-Q 


C = 5rr 


1 + A\ + J~Q 


(27) 


(28) 


(29) 
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II. DYNAMICS OF A SMALL FRAGMENT 
A. Mapping to an effective resisting spring 

The dynamics of a fragment of the filament of length Is ^ (^pas/Gpas) ^ (the strip has lateral extent h) animated by 
the activity of motors may be mapped on to that of a rigid segment working against an effective spring that has both 
passive and active components - here we combine these contributions together in terms of an effective spring constant 
Kg. The passive elasticity of this spring comes from both Gpas and Apas while the active contribution comes from 
attached motors in neighbouring segments. If pn = 0 then the passive part of Kg oc Kpg,s/{isb), while a contribution 
solely due to passive linkers would suggest Kg oc Gpas- Since, active contributions merely rescale we treat this as 
a external parameter. 


B. Degenerate case of free fragment 

When there is no restraining spring, Kg = 0, and equation (11) of the main text is irrelevant. We are left 
with equations (2) and (3) of the main text that describe an untethered, stiff filament forced by motors to move at 
velocity dfU = V. Treating (l^/^m<^on) = ^ as a continuation parameter and restricting ourselves to constant Z, 
we numerically obtain solutions defined by the set N{Z) and Y{Z). Fig. 1 illustrates the dynamical characteristics 
of these solution branches when Ai = ^ and T = exp(f|Y|). The fraction of attached motors, N and Y selects 
the constant speed Z. The stability of the translating solution is controlled by the form of the detachment rate and 
characteristic energies for motors to detach £. The range over which a stable, steadily translating velocity (solid curve) 
can exist is controlled by £. These results are relevant to microtubule-motor assays in which small rigid segments are 
moved about by the motor activity. 



( 35 ) 
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we find three growth rates of the linearized equations of the form —si (real, negative growth rate - the dissipative 
effect), —52 and its complex conjugate, — 52 . The critical frequency at onset uJc = ^c/^on neutral stability 

curve are 


/A/3o + By 

vi+s7Vo; 


(36) 


and 


BYo^T'.No = ^Ai + (1 + + BNo{l + (37) 

Invoking various quasi-static assumptions - thus yielding approximations where the motor displacement is slaved to 
the density field - does not yield qualitative differences. 

The physical mechanism for oscillations is an effective negative spring constant and negative friction co-efficient 
due to the interplay between the motor dynamics and the spring resistance of the rod leading to^ a positive feedback 
in the dynamics of the system. To formally see this we set u = AiU exp {iut) + c.c, T = Yq AiY exp {iut) + c.c and 
N = No-\-AiN exp {iut) and find the find the frequency dependent response of the active force term Fa = pm^m^mYN 
which in frequency space writes as The frequency dependence of y indicates that under certain 

conditions the elastic in-phase - Tl{x) - and/or viscous out of phase - T{x) - parts turn negative suggesting a positive 
feedback thus yielding eventually stable oscillations. For large values of T, the number of attached motors is insufficient 
to supply the energy for oscillations thus yielding an upper bound, for small T too many attached motors increase the 
net elastic resistance and friction coefficients, causing a strong damping of oscillations, these trends are qualitatively 
consistent with related observations on mitotic spindle oscillations even though the governing equation differs due 
the exact averaging procedure [12]. In the stable oscillating state, the active energy input by motor attachment 
must balance the dissipation due to motor friction. Unfortunately, the equation for the emergent frequency is too 
algebraically complicated to be interpreted directly. It is however possible to estimate the effective active viscosity 
and anticipate how this affects emergent frequency. When motors adiabatically couple to the dynamics of the rod, as 
in m for instance, so that Y « ^YF and mean motor extension is slaved to the velocity of the rod and the mean 
time for which the motor stays attached, we find 

,,2 

UJc ^ -=- 

PmkmAl^o{d{i/) I dZ^ 

yielding an estimate of the effective motor friction 

^active ^ Pm^m N,{d{y)/dZ)z^o 

As T increases with everything else held constant, {d{y)/dZ )typically decreases, as does No, resulting in an 
effective increase in ujc. 


(38) 

(39) 


D. Limiting case: = 0 

The ensemble averaged equations for the fraction of attached motors N and the mean motor stretch Y are 

iVi=<n-A'Kn+<ff-^) (40) 


Y, = ^+uj^,^{As-Y) 

Om. 




(41) 


The stability of the stationary base state against small perturbations Z*, T* and A^* is now governed by the evolution 
of the dimensionless set 


Z* = (^N^Y* + YoN*) 


Y* = Z* — Y* 


1-Na 

No 


AT* = -(1 + '^To)N* - ^NoT'oY*. 


(42) 

(43) 

(44) 
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Seeking possible periodic solutions, we again find that the resulting equations admit either decaying solutions or 
oscillatory solutions. The neutral stability curve is given by 

0 = (1 + + BNo{l + 


and the frequency at onset is 


_ _ f^To{iY^:Fo)Y 
V (1 + SiV,) ) 

Again, further simplifications can be made in the various adiabatic limits discussed earlier with motor dynamics 
slaved to the dynamics of the rod. One such physically relevant limit is TJv? ^ 1 which is obtained when the 
detachment rate at the equilibrium point is much larger than the strain-independent attachment rate. An alternate 
approximation which may be treated as a variation of the mean field theory is T ~ As T increases with 

everything else held constant. No decreases, resulting in an effective increase in ujc- 


III. DECAY LENGTH CALCULATIONS 


A. Decay length formula 


Trying out complex solutions of the form U = we get the solvability condition 


(Gpas + Gact^ “ ^pas^ ) + (Gact^) — 0 


Thus A is in general complex and satisfies 

A-i = 


Gi 


pas Gact 


K K 

-f'-pas -f'-pas 




Ga 




pas 


where without loss of generality we have chosen the positive sign for the imaginary part. If 

Gpas , Gact^/ \ 1 Gact \ 

ifi = —— + ——7^(cJ), and ip 2 = — —X(cj), 

-^pas -^pas -^pas 


we get the real and imaginary parts 


Re[XY = ± 




and Im[A^ ] = ± sgn(X) 




(46) 


(47) 


(48) 


(49) 


When the motor aggregates are in rigor^ there is no intrinsic motor dynamics as they are stuck in fixed configurations 
(either attached or detached) and therefore ujon = 0, cjoff = 0 and Vq = 0. In this limit, the motors are bound 
permanently to the filament track and can actually sense the instantaneous displacement of the attachment point 
(where the head is) rather than just the rate of displacement. The number of attached motors N cannot change and 
equals Aq, the preexisting attached fraction. As a result, the Y responds passively to the change in the attachment 
point, and satisfies SmYt = —Uf. Integrating this once yields Sm{Y{t) — Yq) = — Uq). Writing U{t) — Uo = SU 

we get Ew{SU)ss — Pn^n(^G) — kmpmNoi^U) = 0 . Thus the decay length of the strain is given by 

EW \ ^ ^ f -^pas 

kmPmJYo + Pn^N/ VGactA^O + Gpas 

The effective shear modulus due to the links, Geff is obtained by equating the force needed to shear the passive springs 
to the induced extension, (y). That is {NQpmkm + PNkN){dsb){y) = Fg = Gef[{{y)/D)dsb, thus yielding the estimate 
Geff ^ {NoPmkm + PN^n)^^- 





B. Expression for x when As 0 


Recall that Ai = Vokmi^of^Fg) ^ and A 2 = Fs{km^rn) The constant C is then 

^ _ A1A2 f Sm \ 

The (dimensionless) complex susceptibility in this limit is 


where 


Ti = - 



{iuj + + wl^Ai) 


<5n^A°"Uo'A <S^oNo )) 


and 


T 2 


f + ^on + 

\ ^off‘^0^0 / 


Ti. 


C. Expression for y when ^ 0 

Recall that A^ = dj^j The constant C in this limit is evaluated to be 

We find the complex susceptibility is x{^) — (R + 'kF) with 

^ ^ J- 0)2 + + (^ 0^(1 + + ^2 

and 

w ^ ^ g; (1 - UJ <„(1 + ^ 7 b) 

Kff -^ o )2 + + (^ 0^(1 + ^ J -^))2 + ^2 • 
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